First-principles calculations for the adsorption of water molecules on the Cu(lOO) 
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First-principles density-functional theory and supercell models are employed to calculate the ad- 
sorption of water molecules on the Cu(lOO) surface. In agreement with the experimental observa- 
tions, the calculations show that a H2O molecule prefers to bond at a one-fold on-top (Ti) surface 
site with a tilted geometry. At low temperatures, rotational diffusion of the molecular axis of the 
water molecules around the surface normal is predicted to occur at much higher rates than lateral 
diffusion of the molecules. In addition, the calculated binding energy of an adsorbed water molecule 
on the surfaces is significantly smaller than the water sublimation energy, indicating a tendency for 
the formation of water clusters on the Cu(lOO) surface. 

PACS numbers: 68.43.-h, 68.47.-b, 73.20.-r, 82.45.-h 
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I. INTRODUCTION 



The adsorption of water on metal surfaces has been 
extensively studied for many years due to its techno- 
logical importance and wide academic interest iSi^ In 
recent years, there have even been some new experi- 
ments on water enhanced surface reactions of environ- 
mental importancCi^'^ One of the central issues is the 
microscopic understanding of the bonding and orienta- 
tion characteristics of water molecules on the surfaces. 
While experimental and theoretical investigations have 
provided unified information about the bonding struc- 
ture for most water-adsorbed metal surfaces, theoreti- 
cal results are not always in agreement with the experi- 
mental observations. One example is the water-adsorbed 
Cu(lOO) surface&L2i2iiMiii2JS 

Copper is an important material for structural, elec- 
trical, and electronics applications. Corrosion of copper, 
which is a source of significant economic losses and oc- 
curs both in aqueous environments^^'^^ and atmospheric 



air 
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depends intimately on water adsorption. Ad- 
sorption of water on copper is therefore important, both 
from technological and basic scientific points of view. 
The experimental investigations for the adsorption of wa- 
ter on Cu(lOO) have been done with various techniques 
including ultraviolet photoemission spectroscopy (UPS), 
Auger electron spectroscopy (AES), electron-energy- loss 
spectroscopy (EELS), and low-energy electron diffraction 
(LEED)iSiL2ii£ Water has been found to adsorb molecu- 
larly on the surface at low temperatures (< 150 K) and to 
tend to form clusters. Temperature programmed desorp- 
tion (TPD) measurements have shown that the desorp- 
tion of water from the surface occurs at a temperature of 
-160-170 K.^'i" The EELS data of Anderson et al. sug- 
gested that the water molecule is bonded at the one-fold 
on-top (Ti) site with its oxygen end towards the surface 



and its axis tilted away from the surface normal (at an 
estimated angle of about 60°)iSi& 

Several groups have performed theoretical calculations 
for the water-adsorbed Cu (100) surface using cluster 
modelsiii*i2iH While these investigations have provided 
useful insight about the interaction between water and 
the surface, some of the key results are not in agree- 
ment with the experimental data. For example, while 
ah initio density-functional cluster calculations showed 
that the configuration with a tilted water molecule at 
the Ti site on Cu(lOO) would be more stable than that 
with the axis of the water molecule parallel to the surface 
normaljiiii^ they predicted that the Ti site would not be 
the energetically most favorable bonding site. Instead, 
the calculations showed that water would prefer to bond 
at the two- fold bridge (B2) site with its molecular axis 
parallel to the surface normal. -'^'^ Similarly, semi-empirical 
quantum chemical cluster calculations concluded that the 
adsorption of a water molecule at a four-fold hollow site 
{H4) would be energetically most favorable>i2i Both re- 
sults disagree with the experimental data, which suggest 
a tilted geometry of water at the Ti site on Cu(lOO). 

Here we present results of total-energy density- 
functional calculations in which supercell models are used 
to simulate the Cu(lOO) surface. The atomic and elec- 
tronic structures of the adsorbed surface, as well as the 
preferred bonding site and orientation of water molecules, 
are determined. In agreement with the experimental 
data, the most stable adsorption site for water molecules 
chemisorbed on the Cu(lOO) surface is determined to be 
the on-top Ti site, and a tilted geometry of water at the 
on-top site with an angle of approximately 70° away from 
the surface normal is shown to be energetically favorable. 
Moreover, rotational diffusion of the water axis around 
the surface normal is predicted to dominate over lateral 
diffusion at low temperatures. In addition, the energet- 
ics of molecularly adsorbed water on the Cu(lOO) surface 
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suggests a tendency for the formation of water clusters 
on the surface. 

The remainder of this paper is organized as follows. In 
Sec. II we describe the computational method and the 
supercell models that we used. In Sec. Ill we present 
and discuss our theoretical results. Finally, in Sec. IV, 
we summarize the main conclusions obtained from our 
calculations. 



TABLE I: Relaxation of the clean Cu(lOO) surface. Adij is the 
change of the interlayer distance, and do is the corresponding 
distance in the bulk. 





Adia/do (%) 


Ad23/do (%) 


Ad34/do (%) 


This work 


-3.4 


+1.5 


+0.2 


Experiment" 


-2.8 


+1.1 


±0.0 



'Ref. ,3a. 



II. METHOD 

On the (100) surface of an f.c.c. metal, there are three 
different symmetric adsorption sites, known as four-fold 
hollow {H4), two- fold bridge (B2), and one- fold on-top 
(Ti) sitesr^^ We studied the adsorption of water at each 
of the three bonding sites on Cu(lOO), and we also in- 
vestigated the orientation of the water molecule on the 
surface. The Cu(lOO) surface was modeled by repeated 
slabs with seven copper layers separated by a vacuum re- 
gion equivalent to nine copper layers. Each metal layer 
in the supercell contained nine Cu atoms (a 3 x 3 sur- 
face unit cell). Water was adsorbed symmetrically on 
both sides of the slab. All the Cu atoms were initially 
located at their bulk positions, with the equilibrium lat- 
tice constant of the bulk determined by our calculations. 
Calculations with a larger supercell containing nine cop- 
per layers showed that convergence of the binding energy 
of a water molecule on the Cu(lOO) surface for the most 
stable configuration with respect to the thickness of the 
slab was to be within 0.02 eV. 

The calculations were performed within density- 
functional theory, using the pseudopotential method 
and a plane- wave basis setj -^^'^°i^^i^^i^^'^^ The results 
reported in this paper were obtained using the Vi- 
enna ab-initio simulation package (VASP)j22i2ii2^ The 
exchange-correlation effects were treated with the gener- 
alized gradient-corrected exchange-correlation function- 
als (GGA) given by Perdew and Wangi^SsSk We adopted 
the Vanderbilt ultrasoft pseudopotentials supplied by 
Kresse and Hafneri2L2& A plane-wave energy cutoff of 
30 Ry and 10 special k points in the irreducible part 
of the two-dimensional Brillouin zone of the 3x3 sur- 
face cell were used for calculating the H20-adsorbed Cu 
(100) surface. Test calculations with cutoff energies up 
to 35 Ry and with different numbers of special k points 
(6 and 10 points) obtained the binding energy of a wa- 
ter molecule on the Cu(lOO) surface for the most sta- 
ble configuration within 0.01 eV. Optimization of the 
atomic structure was performed for each supercell via a 
conjugate-gradient technique using the total energy and 
the Hellmann-Feynman forces on the atomsii^ All the 
structures were fully relaxed until the forces on all the 
atoms were less than 0.03 eV/A. In this sense, the O-H 
bond length and the H-O-H angle of the H2O molecule 
were optimized. 



III. RESULTS AND DISCUSSIONS 

We first present the calculated properties for bulk cop- 
per and the relaxed clean Cu(lOO) surface. Calculations 
for bulk Cu were conducted with 408 special k points 
and a cutoff energy of 30 Ry. The total energy conver- 
gence with respect to the cutoff energy (25-40 Ry) and 
the number of the k points (220-570) was within a few 
tenths of 1 meV per atom. We obtained lattice constants 
of 3.64 A for bulk Cu, in good agreement with the exper- 
imental value of 3.61 AiSS 

The properties of the clean Cu(lOO) surface were calcu- 
lated using a 1 X 1 surface cell and 66 special k points in 
the surface Brillouin zone. All the seven copper layers ex- 
cept for the central one were relaxed. The surface energy 
of the Cu(lOO) surface was calculated to be 0.59 eV, and 
the results for the surface relaxation are shown in Table 1. 
The surface shows an inward relaxation (approximately 
0.03 A) of the top layer and an outward relaxation (also 
approximately 0.03 A) of the second layer. Our results. 




(a) (b) 




(c) (d) 
Q C» O O OH 

FIG. 1: Schematics of the minimum-energy structures (top 
view) corresponding to the configurations with a water 
molecule at (a) the Ti site (with a tilted geometry), (b) the 
Ti site (without a tilted geometry), (c) the B2 site (with- 
out a tilted geometry), and (d) the H4 site (without a tilted 
geometry) on the Cu(lOO) surface. 
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TABLE II: The binding energies (in unit of eV) of the wa- 
ter molecule adsorbed in the different configurations on the 
Cu(IOO) surface, t/i is the angle between the surface normal 
and the axis of the molecular water. 



Ti {^p = 73°) 


Ti (V' = 0°) 


B2 (^ = 0°) 


H4 ii^ = 0°) 


0.25 


0.17 


0.14 


0.04 



in good agreement with LEED measurementsr^ show re- 
laxation of the Cu(lOO) surface with Adi2/do = —3.4% 
(contraction) and Ad23/c^o = +1.5% (expansion), where 
A(ii2 and Ad23 a-i'e the changes in spacing between the 
top and the second layer and between the second and the 
third layer, and do is the bulk interlayer distance. 

Table II shows the calculated binding energies of a wa- 
ter molecule in the different configurations of the II2O- 
adsorbed Cu(lOO) surface. The most stable configura- 
tion, shown in Fig. 1(a), is that with a tilted water 
molecule at the on-top Ti site (tilted- Ti configuration). 
The water molecule is bonded to the surface with the oxy- 
gen end toward a surface Cu atom and the molecular axis 
tilted away from the surface normal at an angle of 73°. 
The binding energy of a water molecule in this configura- 
tion is calculated to be 0.25 eV. Three metastable config- 
urations (local energy minima) were determined from our 
calculations: a water molecule with its axis parallel to the 
surface normal at the Ti (Fig. 1(b)), B2 (Fig. 1(c)), and 
H4 (Fig. 1(d)) sites, respectively. The binding energies of 
a water molecule in these metastable configurations are 
smaller by 0.08-0.21 eV than that in the tilted- Ti configu- 
ration. We also optimized the structures of the tilted- ,62 
and tilted- configurations (a tilted water molecule at 
the B2 and H4 sites, respectively), and found that nei- 
ther was a local energy minimum. In fact, the tilted 
water molecule in such configurations was found to re- 
lax eventually to a nearby Ti site with a tilted geometry, 
that is, the tilted-Ti configuration (Fig. 1(a)). We there- 
fore conclude that water adsorbed on the Cu(lOO) surface 
prefers the Ti site with its axis tilted away from the sur- 
face normal. This conclusion is in agreement with the 
experimental observations. 

Previous theoretical studies employing cluster models 
also determined the preferred bonding sites of water on 
the Cu(lOO) surface, as mentioned in Sec. I. Ab initio 
DFT calculations using a 13-atom cluster for the Cu(lOO) 
surface showed that water would prefer to bond at the B2 
site with its molecular axis normal to the surfaceii^ Semi- 
empirical quantum chemical cluster calculations with a 
9-atom cluster predicted that the binding energy of a wa- 
ter molecule at the H4 site was significantly larger than at 
the Ti site (by approximately 0.7 eV)fi2, These theoretical 
conclusions are inconsistent with both our ab initio su- 
percell calculations and the experimental data. The main 
problem is that the previous calculations were limited to 
small clusters, containing only up to 13 metal atoms. In 
contrast to semiconductors, metals are often poorly rep- 
resented by small clustersi^'^^i^^i^^i^^i^^i^^'^^i^^i^^ Large 



clusters or extended surface models are usually needed 
to simulate real metal surfaces in order to achieve good 
convergence in calculations. 

The obtained structural parameters of the optimized 
geometry for the tilted- Ti configuration show that the 
water molecule keeps essentially the same shape upon ad- 
sorption on the surface with the equilibrium 0-H bond 
length (0.98 A) and H-O-H angle (105.6°) sfightly larger 
than those of a free molecule: 0.97 A and 105.3°, re- 
spectively, determined from our calculations.^*^ The ex- 
perimental values^ of the 0-H bond length and H-O-H 
bond angle for an isolated water molecule are 0.96 A 
and 104.5°, respectively. After the adsorption of a water 
molecule, the surface Cu atom that is bonded to oxygen 
is observed to shift its position up along the surface nor- 
mal (by 0.08 A) and also to have a slight in-plane shift 
of 0.02 A along the direction against the H2O molecular 
axis. The distance between the oxygen and its bonded 
surface copper atom is 2.28 A. The structural changes re- 
flect the fact that there is a considerable (though weak) 
interaction between the adsorbed water molecule and the 
Cu(lOO) surface. 

It is interesting to note that the variation of the total 
energies for the tilted- Ti configurations with the rotation 
of the water molecule around the surface normal (shown 
in Fig. 2) was determined to be less than 0.02 eV. The ac- 
tivation energy for rotation is thus about 0.08 eV smaller 
than the lateral diffusion barrier for a water molecule 
on the Cu(lOO) surface, which was estimated to be ap- 
proximately 0.1 eV (with T1-B2-T1 taken as the diffusion 



rotation axis 
(surface noraial) 




FIG. 2: Schematics of the rotation of the water molecule 
around the surface normal in the tilted Ti configuration on 
the Cu(lOO) surface. Note that the 0-Cu connecting line is 
slightly off (by approximately 0.03 A) the rotation axis and 
rotates around it. 
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pathway) . While there so far have been no experimental 
efforts to observe such rotational diffusion, these ener- 
getic considerations suggest that it should occur at much 
lower temperatures than lateral diffusion. Lateral diffu- 
sion at elevated temperatures has been experimentally 
observed for single water molecules initially adsorbed at 
low temperatures of about 20 K& The occurrence of the 
rotation of the water molecule around a stiff axis on the 
surface is a consequence of the fact that the bonding of 
water with the Cu(lOO) surface is dominated by the oxy- 
gen lone-pair electronic orbital (see below). 

Reliable theoretical values for the water sublima- 
tion energy have been obtained with ab initio DFT 
calculations 424^ Hamann obtained a value of 0.66 eV 
for the H2O ice sublimation energy;^ while Feibelman;^ 
using the same computation package (VASP) as we em- 
ployed, obtained values of 0.67-0.72 eV for the D2O 
ice sublimation energy. Our calculated binding energy 
(0.25 eV) of H2O on Cu(lOO) is thus significantly smaller 
than the water sublimation energy. Therefore, water 
molecules apparently tend to form ice-like clusters on the 
Cu(lOO) surface. This could occur even at low water cov- 
erages when high temperatures make the water molecules 
diffuse on the surface. Experimental measurements have 
demonstrated that water clusters were indeed formed on 
the Cu(lOO) surfaceAiSiiS In particular, only one peak 
was observed at a temprature of approximately 160 K in 
the TPD spectrum of H2O on Cu(lOO), indistinguishable 
from the ice sublimation peakAiS 

The calculated electronic density of states (DOS) for 
the clean and the H20-adsorbed Cu(lOO) surfaces are 
shown in Fig. 3. The main changes of the DOS upon the 
adsorption of water are found to be in the energy region 
lower than about 4 eV below the Fermi level on the DOS 
curve. Two new peaks (Fig. 3(b)), located at —10.07 
and —6.52 eV relative to the Fermi level, respectively, are 
clearly observed. The curve for the projected density of 
states (projected or local DOS) onto the water molecule 
(Fig. 3(c)) indicates that these two peaks correspond to 
the water I&2 and 3ai states. The highest-energy elec- 
tronic state of water, that is, the I61 molecular orbital, 
is buried in the substrate energy band and hardly ob- 
served in the total DOS curve. However, the local DOS 
curve shows that its peak is located at —4.97 eV rela- 
tive to the Fermi level (Fig. 3(c)). Upon adsorption, the 
Ibi and 3ai orbitals of the water molecule are observed 
to shift down in energy by 0.72 eV and 0.29 eV, respec- 
tively, while the I62 orbital remains almost unchanged in 
energy level. The I61 state is also seen to be more de- 
localized when water is adsorbed on the surface. These 
results suggest that water is bonded to the surface via the 
Ibi state (with some contributions from the 3ai state). 
The Ibi state is primarily of oxygen 2p lone-pair charac- 
ter, and a small charge transfer from water to the copper 
surface is thus mostly from the oxygen 2p state (water is 
known as a Lewis base). Such a charge transfer and the 
resulting bonding between the I61 state and the surface 
would be facilitated when the 2p orbital is perpendicular 




Local DOS 
on H,0 




-12 -9 -6 -3 
Energy (eV) 

FIG. 3: Total density of states (DOS) for the clean (a) and the 
water-adsorbed (b) Cu(lOO) surfaces and local DOS on the 
water molecule (c) for the water-adsorbed Cu(lOO) surface. 
The Fermi level is at eV. The calculated energy levels for the 
molecular orbitals of a free water molecule are also indicated 
in (c). 



to the surface and toward a surface atom. Therefore, a 
water molecule would prefer to bond via oxygen at the 
on-top Ti site on the Cu(lOO) surface. Moreover, the 
molecular axis of water would be tilted away from the 
surface normal since the I61 orbital and the 0-H bonding 
orbitals in water tend to form a quasi-tetrahedral geome- 
try. Thus, the tilted- Ti configuration represents the most 
stable structure for the water-adsorbed Cu(lOO) surface. 
In addition, rotation of the water axis around the surface 
normal has essentially no effects on the bonding between 
the Ibi orbital and the surface, resulting in an very small 
energy cost for such a rotation. These qualitative conclu- 
sions, based on the DOS and electronic structures, are in 
excellent agreement with the quantitative results of the 
energetics discussed above. 



IV. CONCLUSIONS 

The theoretical approach of supercell models combined 
with first-principles total-energy DFT pseudopotential 
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methods has reproduced experimental measurements of 
the preferred adsorption geometry for the water-adsorbed 
Cu(lOO) surface. Very small differences in total ener- 
gies among configurations with different orientations and 
bonding sites of the water molecule have been distin- 
guished, ft is shown that the on-top Ti configuration is 
the most stable structure for the H20-adsorbed Cu(lOO) 
surface, and the water molecule is shown to be tilted with 
its axis approximately 70° away from the surface normal. 
The calculations also predict that rotational diffusion of 
the water axis around the surface normal would occur at 
low temperatures. In addition, water molecules are likely 



to form ice-like clusters on the Cu(lOO) surface at high 
coverages or at high temperatures. 
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